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ABSTRACT 

We present Chandra observations of RXJ1347. 5-1145, the most X-ray luminous cluster 
of galaxies known. We report the discovery of a region of relatively hot, bright X- 
ray emission, located approximately 20 arcsec to the southeast of the main X-ray 
peak, at a position consistent with the region of enhanced Sunyaev-Zeldovich effect 
reported by Komatsu et al. (2001). We suggest that this region contains shocked gas 
resulting from a recent subcluster merger event. Excluding the data for the southeast 
quadrant, the cluster appears relatively relaxed. The X-ray gas temperature rises from 
kT ~ 6 keV within the central 25 h§Q kpc radius to a mean value of ~ 16 keV 
between 0.1 - 0.5 Mpc. The mass profile for the relaxed regions of the cluster, 
determined under the assumption of hydrostatic equilibrium, can be parameterized 
by a Navarro, Frenk & White (1997) model with a scale radius r s ~ 0.4 fc^ Mpc 
and a concentration parameter c ~ 6. The best- fit Chandra mass model is in good 
agreement with independent measurements from weak gravitational lensing studies. 
Strong lensing data for the central regions of the cluster can be also explained by 
the introduction of an additional mass clump centred on the second brightest galaxy. 
We argue that this galaxy is likely to have been the dominant galaxy of the recently 
merged subcluster. 

Key words: galaxies: clusters: individual: RXJ1347. 5-1145 - gravitational lensing - 
X-rays: galaxies - cooling flows - intergalactic medium 
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1 INTRODUCTION 

The launch of the Chandra Observatory (Weisskopf et al. 
2000) in 1999 July has provided the first opportunity for 
detailed, spatially-resolved X-ray spectroscopy of clusters 
of galaxies at moderate to high redshifts. The Advanced 
CCD Imaging Spectrometer (ACIS) on Chandra permits di- 
rect, simultaneous measurements of the X-ray gas temper- 
ature and density profiles in clusters and, via the hydro- 
static assumption, the mass distributions, spanning scales 
from r ~ 10 kpc in cluster cores out to r ~ 1 Mpc at the 
detection limits. 

In this paper, we present the first results from Chandra 
observations of RXJ1347. 5-1145, the most X-ray luminous 
galaxy cluster known (Schindler et al. 1995). This cluster 
has been the subject of several previous X-ray, optical and 
Sunyaev-Zeldovich (SZ) effect studies (e.g. Schindler et al. 
1995, 1997; Fischer & Tyson 1997, Allen 1998; Komatsu et 
al. 1999, 2001; Pointecouteau et al. 1999, 2001). Here, we 
present the first measurements of the X-ray temperature 
structure within the cluster and place tight constraints on 
the total mass distribution. We compare our results with 
those from strong and weak lensing analyses. A comparison 
of the Chandra and SZ results, and a determination of the 



Hubble Constant from the combined data sets, is presented 
elsewhere (Schmidt & Allen 2002, in preparation). 

Except where stated otherwise, the cosmological param- 
eters Ho=50 kms -1 Mpc -1 , Q. = 1 and A = are assumed. 
At the redshift of RXJ1347.5-1145 (z = 0.451), an angular 
scale of 1 arcsec corresponds to physical size of 6.81 kpc in 
this cosmology. 



2 OBSERVATIONS 

The Chandra observations of RXJ1347. 5-1145 were carried 
out using the ACIS on 2000 March 05 and 2000 April 29. The 
target was observed in the back-illuminated S3 detector and 
positioned near the centre of node-1 on CCD 7. The light 
curves for both observations were of high quality with no 
strong background flares. The net good exposure times were 
8.9 and 10.0 ks, respectively, giving a total good exposure 
time of 18.9 ks. The focal plane temperature at the time of 
both observations was -120C. 

We have used the CIAO software and the level-2 events 
files provided by the standard Chandra pipeline processing 
for our analysis. Only those X-ray events with grade classi- 
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Figure 1. (a) The raw 0.3 — 7.0 koV Chandra image of RXJ1347. 5-1145. The pixel size is 4 detector pixels (1.97 arcsec). (b) Contour 
plot of the same region, adaptively smoothed using the code of Ebeling et al. (2002) with a threshold value of 3<r. The contours have 
equal logarithmic spacing. 



fications of 0,2,3,4 and 6 were included in our final cleaned 
data sets. 



3 X-RAY IMAGING 
3.1 X-ray morphology 

The raw 0.3 — 7.0 keV image of the central 6x6 arcmin 2 
region of the cluster, from the combined 18.9 ks data set, 
is shown in Fig. |l|(a). The pixel size is 1.97 x 1.97 arcsec 2 , 
corresponding to 4 x 4 raw detector pixels. Fig. [j](b) shows 
an adaptively smoothed contour plot of the same data, us- 
ing the smoothing algorithm of Ebeling, White & Rangara- 
jan (2002). The images reveal a number of notable fea- 
tures. Firstly, we see a sharp central surface brightness peak 
at a position 13h47m30.63 -Ild45m09.3s (J2000.), in good 
agreement with the optical centroid for the dominant clus- 
ter galaxy, 13h47m30.5 -Ild45m09s (J2000.; Schindler et al. 
1995). Secondly, we identify a region of enhanced emission 
~ 20 arcsec to the southeast of the X-ray peak (see below) . 
Thirdly, on large scales (r ~ 80 arcsec) we detect an exten- 
sion of the X-ray isophotes to the south. 

Fig. |^ shows the data for the central regions of the clus- 
ter on a finer 0.492 x 0.492 arcsec 2 scale (corresponding to 
lxl raw detector pixels). From this figure it appears that 
the region of enhanced emission to the southeast of the X- 
ray peak has a roughly circular morphology with a flattened 
northwest edge. Beyond this edge, at a radius r ~ 10 arcsec, 
the X-ray subclump appears to be separated from the main 
cluster core by a valley of reduced emission. The spectral 
analysis in Section 4.2 shows that the X-ray subclump has 
a hotter temperature than the gas at the same radius in 
other directions from the X-ray peak and probably contains 
shocked gas resulting from a recent subcluster merger event. 



Komatsu et al. (2001) report the detection of a region 
of enhanced SZ effect at a position coincident with the X- 
ray subclump. The enhanced SZ effect is consistent with 
the higher temperature and density for this region measured 
from the Chandra data. 

No strong X-ray point sources are detected in the Chan- 
dra S3 field. (The brightest point source has a flux of 
3.1 x 10~ 14 erg cm -2 s _1 in the 0.5 - 7.0 keV band; Gandhi 
et al. , in preparation.) 



3.2 Surface brightness profiles 

The azimuthally-averaged, 0.3 — 7.0 keV X-ray surface 
brightness profile for RXJ1347. 5-1145 for position angles of 
180 — 90 degrees (i.e. excluding the southeast quadrant) and 
for the disturbed, southeast quadrant are shown in Figs. 
la,b), respectively. The profiles have been flat-fielded and 
background subtracted using an on-chip region free from 
cluster emission. The bin-size is 2 detector pixels (0.984 arc- 
sec). 

We see that the data for position angles of 180 — 90 
degrees appear regular and, for r < 500 kpc, can be ap- 
proximated (x 2 = 102 for 72 degrees of freedom) by a f3- 
model (e.g. Jones & Forman 1984) of the form 5x(r) = 
S(0) [l + (r/r c ) 2 ] 1/2-3/3 with a core radius r c = 29.2 ± 0.7 
kpc and a slope parameter /3 = 0.535 ± 0.003 (la errors). 
Between radii of 0.1 — 1.0 Mpc we find evidence for a steep- 
ening of the surface brightness profile with increasing radius: 
for 0.1 < r < 1.0 Mpc the data can be described (\ 2 = 135 
for 131 degrees of freedom) by a simple broken power-law 
model, with a break at a radius of 487^ 2 g kpc, and slopes 
in the regions internal and external to the break radius of 
—2.25 ± 0.03 and —3.66 ± 0.18, respectively. A similar re- 
sult was obtained for Abell 2390 by Allen, Ettori & Fabian 
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Figure 2. A raw 0.3 — 7.0 kcV image of the central regions of 
RXJ1347.5-1145 on a finer spatial scale (pixel size 0.492 X 0.492 
arcsec 2 , equivalent to 1 X 1 raw detector pixels). 



Table 1. The results from the analysis of the annular spectra 
(covering position angles 180 — 90 degrees). Temperatures (fcT) 
are in keV and metallicities (Z) in solar units. The absorbing 
column density has been fixed at the nominal Galactic value of 
4.85 xlO 20 atom cm 2 . The total \ 2 values and number of degrees 
of freedom (DOF) in the fits are listed in column 4. Error bars 
are the la (Ax 2 = 1.0) confidence limits on a single interesting 
parameter. 
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(2001b), although in that case the slope interior to the break 
radius was flatter. 

The data for the southeast quadrant containing the 
X-ray subclump reveal a clear excess of emission relative 
to other directions. Enhanced emission is detected between 
radii of ~ 70 — 300 kpc, with a maximum enhancement of a 
factor ~ 3 at a radius r ~ 180 kpc. 



4 SPATIALLY-RESOLVED SPECTROSCOPY 
4.1 Method 

For our spectral analysis, we have separated the southeast 
quadrant containing the X-ray subclump from the rest of the 
cluster. The data for position angles of 180 — 90 (the relaxed 
part of the cluster) were divided into the annular regions de- 
tailed in Table l.n A single spectrum was extracted for each 
region in 1024 Pulse Invariant (PI) channels. The spectra 
were re-grouped to contain a minimum of 20 counts per PI 
channel thereby allowing x 2 statistics to be used. 

Background spectra, appropriate for the regions stud- 
ied, were extracted from the ACIS-S3 blank-field data sets 
provided by Maxim Markevitch (which are available from 
the Chandra X-ray Center). Separate photon- weighted re- 
sponse matrices and effective area files were constructed for 
each region using the available calibration and response files 
appropriate for the -120C focal plane temperature.^ 

* For the southeast quadrant, we have extracted a single spec- 
trum covering radii 60—195 kpc (18-58 pixels), where the emission 
enhancement in Fig. ^ is most obvious. 

t For each 32 X 32 pixel 2 sub-region of the S3 chip, a spectral 
response (.rmf) and an auxiliary response (.arf) matrix were cre- 
ated using the CIAO tools rakrmf and mkarf, respectively. For 



The data from both Chandra observations were mod- 
elled simultaneously using the XSPEC code (version 11.01; 
Arnaud 1996). We have limited our analysis to the 0.5 — 7.0 
keV energy band, over which the calibration of the back- 
illuminated CCD detectors is currently best understood. 

The spectra have been modelled using the MEKAL 
plasma emission code of Kaastra & Mewe (1993; incor- 
porating the Fe L calculations of Liedhal, Osterheld & 
Goldstein 1995) and the photoelectric absorption models of 
Balucinska-Church & McCammon (1992). We have fitted 
each annular spectrum using a simple, single-temperature 
model with the absorbing column density fixed at the nom- 
inal Galactic value (A+i = 4.85 x 10 20 atom cm~ 2 ; Dickey & 
Lockman 1990). The free parameters in this model were the 
temperature (fcT), metallicity (Z; measured relative to the 
solar photospheric values of Anders & Grevesse 1989, with 
the various elements assumed to be present in their solar 
ratios) and emission measure. We note that including the 
absorbing column density as an additional free parameter 
did not result in significant improvements in the goodness 
of fit. 

4.2 Results 

The best-fit parameter values and la (A^ 2 = 1.0) confi- 
dence limits determined from the fits to the annular spec- 
tra (excluding the southeast quadrant) in the 0.5 — 7.0 keV 
band are summarized in Table 1. The projected temperature 

each of the annular regions studied, the number of source counts 
in each 32 X 32 pixel 2 sub-region was determined. The individ- 
ual .rmf and .arf files were then combined using the FTOOLS 
programs addrmf and addarf (using a script provided by Rod- 
erick Johnstone) to form a counts-weighted spectral response and 
auxiliary response matrix appropriate for each annulus. 
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Figure 3. Background-subtracted, flat-fielded, azimuthally-averaged radial surface brightness profiles for RXJ1347. 5-1145 in the 0.3 — 7.0 
keV band. The bin-size is 0.984 arcsec (6.70 kpc). (a) The results for position angles of 180 — 90 degrees. The profile appears regular 
and relaxed (b) The results for position angles of 90 — 180 degrees (the southeast quadrant). The dotted curve shows the results for the 
180 — 90 degree region. Note the clear excess of emission in the southeast quadrant between radii of 10 — 45 arcsec (70 — 300 kpc). 



profile determined with this model is shown in Fig. ^. The 
temperature rises from a mean value of 7.0l 'g keV within 
r = 24 kpc to kT = 15.8±\-% keV over the 0.1 - 0.5 Mpc 
region. 

A fit with the same model to the data for the southeast 
quadrant between radii of 60 — 195 kpc (where the X-ray 
enhancement is most obvious in the image) gives a best-fit 
temperature kT — 18.0^2 3 keV. (A similar result, kT = 
18. keV, is obtained using a circular region of radius 
10 arcsec centred on the enhancement.) In other directions, 
the mean temperature between radii of 60 — 195 kpc is kT = 
12.7 ± 1 keV. The higher temperature and enhanced X-ray 
surface brightness in the southeast quadrant are consistent 
with the X-ray gas in that region having undergone shock 
compression (Section 8). 

We find marginal evidence for a metallicity gradient in 
the cluster, with a mean metallicity for the central 100 kpc 
radius of Z = 0.55±0.09 solar, as opposed to Z = 0.34±0.12 
solar for the 100 — 740 kpc region. 
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4.3 Spectral deprojection analysis 

The results discussed above are based on the analysis of pro- 
jected spectra. We have also carried out a simple deprojec- 
tion analysis of the Chandra spectral data using the method 
described by Allen et al. (2001b). 

For this analysis we have used the same annular regions 
(covering position angles of 180 — 90 degrees) and have as- 
sumed that the emission from each spherical shell (the shells 
are defined by the same inner and outer radii as the an- 
nular regions) is isothermal and absorbed by the Galactic 
column density. The fit to the outermost annulus is used to 
determine the temperature and emission measure in the out- 
ermost spherical shell. The contribution from that shell to 
each inner annulus is then determined by purely geometric 
factors (e.g. Kriss, Cioffi & Canizares 1983). The fit to the 



Figure 4. The projected X-ray gas temperature profile (and 1<t 
errors) measured from the Chandra data in the 0.5 — 7.0 keV 
energy band. 



second annulus inward is used to determine the parameters 
for the second spherical shell, and so forth, working inwards. 

The data for all eight annular spectra were fitted si- 
multaneously in order to determine the parameter values 
and confidence limits. The metallicity was linked to take the 
same value at all radii. The temperature profile determined 
with the spectral deprojection method is shown in Fig. n. 



© 0000 RAS, MNRAS 000, 000-000 



Chandra observations of RXJ1347. 5-1145 5 



o 

CV2 




10 100 
Radius (kpc) 

Figure 5. The predicted deprojected temperature profile (grey 
curve) determined from 100 Monte-Carlo simulations using the 
best-fitting NFW mass model (with r s = 0.40 Mpc, c = 5.87 
and a = 1450 km s — 1 ; Section 5). The predicted profile has been 
binned to the same spatial resolution as the spectral deprojcction 
results (solid points; Section 4.3), which have been overlaid. The 
agreement between the deprojected spectral results and the best- 
fit NFW mass model predictions (reduced xt = 0.87 for v = 6 
degrees of freedom) indicates that the NFW mass model provides 
a good description of the spatially-resolved Chandra spectra. 

4.4 Comparison with previous work 

The mean emission-weighted temperature and metallicity 
for RXJ1347.5-1145 of kT = 12.2 ± 0.6 keV and Z = 
0.42 ± 0.07 solar, respectively, measured in the 0.0 — 0.74 
Mpc range over the full 360 degrees, are in good agree- 
ment with the values reported by Allen & Fabian (1998; 
kT = 12.5±g;| keV, Z = 0.38±J^J solar) using ASCA data 
and a similar spectral model. Excluding the data for the 
southeast quadrant, we obtain only small changes in these 
results: kT = 12.0 ± 0.6 keV and Z = 0.41 ± 0.07 solar. 
The Chandra results are also in good agreement with those 
reported by Schindler et al. (1997) using the same ASCA 
data and a 3 arcmin radius aperture (kT = 11.8ICi'o keV, 
Z = 0.39^0 solar). The emission- weighted temperature 
measured with Chandra is slightly cooler than the value of 
14.5ti g keV measured with BeppoSAX by Ettori, Allen & 
Fabian (2001), as can be expected given the observed tem- 
perature profile (Figs. Q^) and the different response char- 
acteristics of the instruments. (The data for the southeast 
quadrant are not excluded from the BeppoSAX or ASCA 
studies.) 

The 'ambient' emission- weighted temperature measured 
with Chandra in the 0.1 — 0.5 Mpc range (i.e. excluding the 
data for the southeast quadrant, the central, cool region, and 
the regions at large radii where the temperature may start 
to decline again) of kT — 15.8 jz^"® keV is in good agreement 
with the value of 15.9 6 keV estimated from BeppoSAX 
data by Ettori et al. (2001) using a spectral model which 
included a constant-pressure cooling flow component. The 
Chandra result is also consistent with the ambient temper- 
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Figure 6. Contour plot showing the 68.3% (la), 95.4% (2a) and 
99.73% (3a) confidence contours on the scale radius, r s , and the 
effective velocity dispersion, a, for the NFW mass models. Con- 
tours of constant concentration parameter, c, are marked with 
dashed lines. The best-fit model is marked by a plus sign. 

ature of kT = 26lf 2 keV estimated from ASCA data us- 
ing a similar constant-pressure cooling flow model (Allen & 
Fabian 1998), and kT = I8.6I3 J keV from a re-analysis of 
the ASCA data by the same authors using the more appro- 
priate 'isothermal' cooling flow models of Nulsen (1998). 

The emission weighted ambient temperature deter- 
mined from the Chandra data in the 0.1 — 0.5 Mpc region is 
in good agreement with the mean gas mass-weighted tem- 
perature of I6.1I2 7 keV measured within r2soo (r = 0.72 
Mpc) by Allen, Schmidt & Fabian (2001c). We note that 
this temperature is also consistent with the predicted value 
of kT = 16.8 keV, using the cooling-flow corrected fcTx/^Boi 
relation of Allen & Fabian (1998; we assume a bolometric 
luminosity of ~ 2.2 x 10 46 ergs -1 as measured by ASCA). 



5 MEASUREMENT OF THE CLUSTER MASS 
PROFILE 

5.1 Method 

The observed X-ray surface brightness profile (Fig. |^a) and 
deprojected, spectrally-determined temperature profile (Fig. 
|E]) may together be used to determine the X-ray gas mass 
and total mass profiles in the cluster. For this analysis, we 
have used an enhanced version of the image deprojection 
code described by White, Jones & Forman (1997), and have 
followed the methods outlined by Allen et al. (2001b) and 
Schmidt et al. (2001). 

A variety of simple parameterizations for the cluster 
mass distribution were examined to establish which could 
provide an adequate description of the Chandra data (see 
below). The best-fit parameter values and confidence limits 
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were determined by examining parameter grids and evaluat- 
ing x 2 f° r each set of parameters.^ Spherical symmetry and 
hydrostatic equilibrium are assumed throughout. 



5.2 NFW mass models 

We find that a good fit (Xmin = 5-2 for 6 degrees of freedom, 
hereafter DOF) to the Chandra data for RXJ1347.5-1145 
can be obtained using a Navarro, Frenk & White (1997, 
hereafter NFW) model: 



p(r) = 



Pcrit(z)5c 



(r/r s ) (l+r/r s y 



(1) 



where p(r) is the mass density, p cr it{z) = 3H(z) 2 /8ttG is the 
critical density for closure and 



5c 



200 



3 [m(l + c) - c/(l + c)] 



(2) 



The best-fit scale radius, r B = 0.40_ 012 Mpc and the con- 
centration parameter, c = 5.87tJ'44 (68 per cent confi- 
dence limits). The normalization of the mass profile may 
also be expressed in terms of an effective velocity disper- 
sion, a = \/^0H (z)r s c = 1450l2oo kms -1 (with r B in units 
of Mpc). The deprojected X-ray gas temperature profile pre- 
dicted by the best-fitting NFW mass model (given the ob- 
served surface brightness profile) is shown overlaid on the 
observed spectral results in Fig. q. 

Fig. ^ shows a contour plot of the \ 2 values obtained 
for the full range of NFW models studied. The minimum 
X 2 value is marked with a cross. Contours have been drawn 
at intervals of Ax 2 = 2.30, 6.17 and 11.8, corresponding to 
formal confidence limits of 68.3 (If), 95.4 (2a) and 99.73 
(3<r) per cent for two interesting parameters. Fig. ^ shows 
the best-fitting NFW mass profile, which has a virial radius 
?"200 = cr B — 2.351q 33 Mpc and an integrated mass within 
the virial radius, M 20 o = 2.29+J;^ x 10 15 M . 



5.3 Other parameterized mass models 

We have also examined a variety of other popular parameter- 
ized mass models. Firstly, the model of Moore et al. (1998): 

i<5 c „ 



p(r) 



where 



Pcrit(z)' 



(r/ry) 1 - 5 (1 + r/r s y 



(3) 



t The observed surface brightness profile (Fig. |i]a.) and a partic- 
ular parameterized mass model are together used to predict the 
temperature profile of the X-ray gas. (We use the median tem- 
perature profile determined from 100 Monte-Carlo simulations in 
which the surface brightness profile is perturbed according to the 
statistical uncertainties at each radius. The outermost pressure is 
fixed using an iterative technique which ensures a smooth pres- 
sure gradient in these regions.) The predicted temperature profile 
is rebinned to the same binning as the spectral results and the 
X 2 difference between the observed and predicted, deprojected 
temperature profiles is calculated. The parameters for the mass 
model are stepped through a regular grid of values in the r B -cr 
plane (see text) to determine the best-fit values and 68 per cent 
confidence limits. The gas mass profile is determined to high pre- 
cision at each grid point using the observed surface brightness 
profile and model temperature profile. 
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Figure 7. The integrated (3-dimensional) mass profile and la 
errors for RXJ1347. 5-1145 determined from the Chandra data 
using the NFW parameterization (section 5.2). The grey curve 
shows the X-ray gas mass profile. (The best-fit X-ray gas mass 
profile and la errors are plotted, although these curves are indis- 
tinguishable in the figure.) 



100 



ln(l + c m 5 ) ' 



(4) 



For this model we obtain x min = 8.6 (6 DOF), with 
best-fit values for the scale radius r B — 5.0 Mpc (the 
maximum allowed value in our grid), concentration pa- 
rameter c m = 0.69 and effective velocity dispersion a = 
2125 kms -1 (where a = ■s/50Ho(z)r e c m )- 

Secondly, we examined a simple, non-singular isother- 
mal sphere model: 



p(r) = 



1 



2ttG 



+ ri 



(5) 



for which we obtain Xmin = 4.9 (6 DOF) with r c = 45 ± 10 
kpc and o"i so = 1590 ± 150 kms -1 . 

We thus find that all of the two-parameter mass mod- 
els described above provide acceptable descriptions of the 
Chandra data for RXJ1347. 5-1145. We note, however, that 
a singular isothermal sphere [p(r) oc r~ 2 ] does not provide 



an acceptable fit (x 



2 
min 



33.1 for 7 DOF). 



6 COMPARISON WITH GRAVITATIONAL 
LENSING RESULTS 

6.1 Weak lensing 

Fischer & Tyson (1997) present a detailed weak lensing 
study of RXJ1347. 5-1145. They measure a mass profile that 
can be parameterized by an NFW model with r2oo = 2.4 
Mpc, or a singular isothermal sphere with a = 1500 kms -1 , 
in good agreement with the Chandra results. 

The azimuthally-averaged, projected surface mass den- 
sity profile for RXJ1347. 5-1145 determined by Fischer & 
Tyson (1997) is shown in Fig. ^. We have excluded the in- 
nermost point from their data, which extends into the strong 
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Figure 8. A comparison of the projected surface mass density contrast determined from the Chandra X-ray data (Section 5) with the 
weak lensing results of Fischer &; Tyson (1997; solid triangles) and the strong lensing result from Section 6.2 (grey circle). The best-fit 
NFW X-ray mass model for the relaxed regions of the cluster and la confidence limits (the maximum and minimum values at each radius 
for all NFW models within the 68 per cent confidence contour shown in Fig. ||) are shown as the dotted and full curves, respectively. The 
strong lensing point is marked with a grey circle. Note that the strong lensing point should lie above the X-ray results, which exclude 
the regions of the cluster affected by the second mass clump. We adopt r max = 2.72 Mpc as in Fischer & Tyson (1997). 



lensing regime (see below; r max = 2.72 Mpc in our adopted 
cosmology) . The same quantity, calculated for all NFW mass 
models within the 68 per cent confidence contour determined 
from the Chandra X-ray data (Fig. ^) has been overlaid. We 
see that the Chandra and weak-lensing results exhibit good 
agreement on all spatial scales studied. 

The agreement between the independent lensing and 
X-ray mass measurements confirms the validity of the hy- 
drostatic assumption used in the X-ray analysis (having ex- 
cluded the southeast quadrant) and suggests that the mass 
profile in the cluster has been robustly determined. 



6.2 Strong lensing 

Schindler et al. (1995) report the discovery of two bright 
gravitationally-lensed arcs in RXJ1347. 5-1145 using ground 
based optical imaging. Sahu et al. (1998) present Hubble 
Space Telescope Imaging Spectrograph (STIS) observations 
of the cluster which reveal a number of additional strongly- 



lensed features. The STIS image for RXJ1347.5-1145 is 
shown in the left panel of Fig. M. Sahu et al. (1998) also 
present ground-based spectroscopy of the bright northern 
arc in the cluster (arc 1) at a radius of 34.9 arcsec (240 kpc) 
from the dominant cluster galaxy, for which they measure a 
redshift zi = 0.806. 

Following the methods outlined in Schmidt et al. (2001), 
we have examined the constraints that the observed strong 
lensing configuration can place on the mass distribution in 
the cluster core. We first examined a simple, spherically- 
symmetric mass model, centred on the dominant cluster 
galaxy, with the best-fit parameter values determined from 
the Chandra data for the relaxed regions of the cluster 
(Sect. 5.2). We find that such a mass model does not pro- 
vide a sufficient central mass density to explain the strong 
lensing data. 

We next examined refined lensing models in which a 
second mass component, centered on the second brightest 
galaxy was introduced, together with appropriate elliptici- 
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Figure 9. The (a) observed and (b) predicted gravitational arc geometry in RXJ1347. 5-1145. The image in the left panel has been 
compiled from archival HST STIS data that were originally presented by Sahu et al. (1998). The arrow points to an additional arc-like 
feature mentioned in the text. The right panel shows the arc geometry predicted by our two-component mass model. The arc numbering 
has been taken from Sahu et al. (1998). The arcs have been simulated using circular Gaussian surface brightness distributions. The 10%, 
30% and 50% brightness contours are plotted. The true source positions are denoted by crosses, with details in Table ^. The central 
positions for the two mass components are marked with plus signs. Also shown are the critical curves (dash-dotted) and caustic lines 
(solid; e.g., Schneider, Ehlers & Falco 1992) for a source at a redshift of z = 0.97. 



ties for the matter distributions. (The introduction of such 
a mass clump is well-motivated by the Chandra data for 
the southeast quadrant of the cluster. The second bright- 
est galaxy is likely to have been the dominant galaxy of the 
recently- merged subcluster; Section 8). Both mass compo- 
nents have been modelled by elliptical NFW models (see 
Schmidt et al. 2001 for definitions). The parameters for the 
main mass component were fixed to the values determined 
from the Chandra data for the relaxed regions of the cluster 
in Sect. i.e. scale radius r B = 0.4 Mpc, velocity disper- 
sion a — 1450 kms" 1 and concentration parameter c = 5.87. 
The axis ratio and position angle were fixed to the values 
measured for the dominant galaxy from the STIS data: axis 
ratio q = 0.76 and position angle 9 = 13.1°. The centre 
of the main mass component was fixed to the centre of the 
dominant galaxy. 

The parameters for the second mass component have 
been determined by requiring that the overall potential 
is able to produce the northern arc (arc 1) at the mea- 
sured redshift, z\ = 0.806. For an assumed scale radius 
r Sl ciump = 0.25 Mpc, we obtain a best-fitting velocity dis- 
persion (Tdump = 815 km s" 1 and a concentration param- 
eter c c i um p = 5.3. Such parameters are typical for mid- 
temperature galaxy clusters (kT ~ 5 keV. e.g. Allen et 
al. 2001c). The position, ellipticity and position angle of 
the second mass component were fixed to the values deter- 
mined for the second brightest galaxy from the STIS image: 
gdump = 0.79 and # c i ump = 28.1°. 

The redshift measurement for the northern arc con- 
strains the projected mass within the arc radius. It is possi- 



Table 2. The redshifts, sizes (in arcsec) and positions (offsets in 
arcsec with respect to the dominant cluster galaxy in the STIS 
reference frame) for the sources in the lens model shown in Fig. Pi 



arc 


redshift 


size 


position 






(FWHM) 


A x 


Ay 


1 


0.806 


0.8 


2.77 


10.03 


2 


0.75 


0.6 


5.71 


9.44 


3 


0.97 


0.6 


1.00 


-6.23 


1 


0.97 


0.5 


-7.24 


-2.66 


5 


0.97 


0.(5 


-13.44 


-1.16 



ble to vary the scale radius for the second mass component 
and to adjust the concentration parameter accordingly, so 
long as the projected mass within the arc radius remains 
constant. The absence of multiple arc images within the 
cluster core suggests r s ,clump > 0.1 Mpc. 

Having adjusted the parameters for the second mass 
component to produce the northern arc, we then examined 
whether the resulting mass model can explain the other arcs 
observed by Sahu et al. (1998). Unfortunately, there are no 
measured redshifts for the other arcs. (Sahu et al. 1998 note 
a possible lower limit for arc 4 of z > 1 based on the non- 
detection of [On] in their ground-based spectra, although 
they caution that the arc may simply be a galaxy with- 
out a strong emission- line feature). The details of our model 
predictions are shown in the right panel of Fig. ^| The re- 
sults on the predicted arc redshifts and the assumed sizes 
of the sources are listed in Table H. The predicted redshifts 
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are cosmology dependent and would be modified slightly 
if, for example, an accelerating universe with a different 
Hubble Constant were assumed. For the mass model used 
here, the predicted redshifts should be accurate to ~ ±0.05, 
with the exception of arc 4 which is determined more pre- 
cisely. We have modelled the lensed galaxies using circu- 
lar Gaussian brightness distributions with the full-width- 
at half-maximum (FWHM) values listed in Table ^. The 
FWHMs for the arcs were adjusted so that their 30 per cent 
brightness contours provide a good fit to the STIS data. 
The effects of using different source sizes is illustrated by 
the different brightness contours in Fig. ^[ 

Although we include arc 2 in our modelling, we caution 
that this source appears quite red in the colour image pub- 
lished by Fischer & Tyson (1997) and that it is possible that 
the source is simply an edge-on spiral galaxy. The predicted 
redshift for arc 2 is lower than zi, since otherwise the lens 
model would predict a much stronger distortion. Arcs 3, 4 
and 5 are consistent with being at the same redshift. The 
projected separation of the sources producing arcs 4 and 5 
is 43.4 h^g kpc. Although arcs 4 and 5 have a similar colour 
in the data of Fischer & Tyson (1997), it is not possible to 
explain them as images of the same source using our mass 
model. In exploring different redshifts for arc 3, we found 
that for redshifts higher than 2:3 ~ 1 it is possible to produce 
a counter image on the other side of the central galaxy. This 
alerted us to the presence of a blue object (identified from 
the colour image of Fischer & Tyson 1997) to the south-east 
of the merging subclump (see the arrow in the left panel of 
Fig. ^). However, both the appearance in the STIS image 
and our lensing analysis suggest that these objects are not 
images of the same source. 

We note that circular NFW models for one or both of 
the main mass components cannot explain the observed arc 
shapes and orientations well. In particular, the orientation of 
arc 4 reflects the orientation of the critical curve (Fig. ^) in 
its vicinity. A small ~ 5° difference between the orientation 
of the observed and predicted arc 4 remains in our best-fit 
model, which points to a small imprecision of the lensing 
model. 

The grey circle in Fig. ^ shows the azimuthally- 
averaged, surface mass density contrast at r = 240 kpc 
(centred on the main mass component) for the best-fit strong 
lensing model. Note that the strong lensing results should lie 
above the X-ray results at the same radii, since the Chandra 
analysis excludes the regions of the cluster affected by the 
second mass subclump. The surface mass density contrast 
rises sharply between the innermost weak lensing point and 
the strong lensing result, suggesting that the second mass 
component is centrally concentrated and does not make a 
major contribution to the total mass beyond the strong lens- 
ing regime. 

In conclusion we find that a simple two-component mass 
model, with ellipticities and orientations for the mass com- 
ponents matching those of the dominant cluster galaxies, 
provides a reasonable description for the overall arc geome- 
try in RXJ1347. 5-1145. Such a model is consistent with the 
Chandra X-ray and weak lensing results. 



7 THE PROPERTIES OF THE CLUSTER GAS 



7.1 Electron density and cooling time profiles 

The results on the electron density and cooling time as a 
function of radius, determined from the image deprojection 
analysis using the best-fit NFW mass model, are shown 
in Fig. |Io| Within the central 500 kpc radius, the elec- 
tron density profile can be parameterized (\ 2 = 43.2 for 
34 degrees of freedom) by a /3-profile with a core radius, 
r c = 27.8 ± 1.6 kpc, (3 = 0.521 ±0.012 and a central density, 
n e (0) = 2.3 ± 0.1 x 10~ 2 cm~ 3 (la errors). We measure a 
central cooling time (i.e. the mean cooling time within the 
central 1.97 arcsec (~ 13.4 kpc) bin of t coo i = 4.1±gJ x 10 8 
yr (uncertainties are the 10 and 90 percentile values from 
100 Monte Carlo simulations). 



7.2 The X-ray gas mass fraction 

The X-ray gas-to-total-mass ratio as a function of radius, 

fsas(r), determined from the Chandra data is shown in Fig. 

We find that the best-fit / gas value rises rapidly with 

increasing radius within the central r < 40 kpc radius and 

then remains approximately constant, or rises slowly, out to 

the limits of the data at r = 0.74 Mpc where we measure 
f - n 1 41 +0.035 

Following the usual arguments, which assume that the 
properties of clusters provide a fair sample of those of the 
Universe as a whole (e.g. White & Frenk 1991; White et al. 
1993; White & Fabian 1995; Evrard 1997; Ettori & Fabian 
1999), we may use our results on the X-ray gas mass frac- 
tion to estimate the total matter density in the Universe, 
fl m . Assuming that the luminous baryonic mass in galax- 
ies in RXJ1347-1145 is 0.134 h% 5 of the X-ray gas mass 
(e.g. White et al. 1993; Fukugita, Hogan & Peebles 1998) 
and neglecting other possible sources of baryonic dark mat- 
ter in the cluster, one can show that fl m = f2b/1.134/ gas , 
where fib is mean baryon density in the Universe. For 
fib/150 = 0.0820 ± 0.0072 (O'Meara et al. 2001) we ob- 
tain fi m = 0.51 ± 0.12. For a ACDM cosmology with Ho = 
70kms _1 Mpc -1 , our results for RXJ1347.5-1145 would im- 
ply fi m = 0.33 ± 0.08. 



7.3 Cooling flow models 

Under the assumption that the X-ray gas in the core of 
RXJ1347. 5-1145 is in a steady-state cooling flow, we can 
parameterize the luminosity profile in terms of an equivalent 
mass deposition rate (e.g. White et al. 1997). The results 
of this calculation are shown in Fig. [l^. We see that the 
integrated mass deposition rate (MDR) rises with increasing 
radius within the central 100 kpc, to a maximum value of 
~ 1300 M Q yr" 1 . 

We have examined the spectrum for the central 100 
kpc (radius) region using a model appropriate for a cooling 
flow with distributed mass deposition. (We use the model 
of Nulsen 1998 and assume that the integrated MDR within 
radius r, M oc r.) The normalization of the cooling-flow 
component is parameterized in terms of a mass deposition 
rate, M. The cooling flow was also allowed to be absorbed 
by an intrinsic column density, ANn, of cold gas (with solar 
metallicity) located at the redshift of the cluster. Both M 
and ANh were free parameters in the fits. Accounting for 
projection effects, we find that the 0.5 — 7.0 keV spectrum 
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Figure 10. The results on (a) the electron density and (b) the cooling time, determined from the X-ray image dcprojection analysis 
using the best-fit NFW mass model. Error bars are the la errors determined from 100 Monte Carlo simulations. A Galactic column 
density of 4.85 X 10 20 atom cm -2 and a metallicity of 0.4 solar are assumed. 



C\2 




Radius (kpc) 

Figure 11. The ratio of the X-ray gas mass to total gravitating 
mass as a function of radius. The three curves show the best-fit 
value (dotted curve) and la confidence limits (solid curves). At 
r = 0.74 Mpc we measure / gas = 0.141 j;2;™. 
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100 200 

Radius (kpc) 



Figure 12. The mass deposition rate (MDR) determined from 
the image deprojection analysis under the assumption that the 
central regions of the cluster contain a steady-state, inhomoge- 
neous cooling flow. Error bars are the 10 and 90 percentile values 
from 100 Monte Carlo simulations. The dotted line shows the 
best-fitting broken power-law model described in Section 7.3. 



for the central 100 kpc radius region is consistent with a 
cooling flow of < 600 Mq yr -1 , intrinsically absorbed by 
a column density of 1 — 4 x 10 21 atom cm -2 , with a mean 
mass-weighted flow temperature of ~ 10 keV. 

Fitting the mass deposition profile in Fig. jl2] with a 
broken power-law model, we determine a break radius of 
31 ± 2 kpc, and slopes internal and external to the break 
radius of 2.7 ± 0.3 and 0.5 ± 0.1, respectively (la errors). 
Allen et al. (2001a) argue that such breaks may mark the 



outer edges of the present-day cooling flows in clusters. The 
predicted MDR for the central r ~ 20 kpc radius, where a 
steady cooling flow with gas cooling to zero degrees might 
be expected to occur, is ~ 200 M© yr -1 . The cooling time 
of the X-ray gas at the break radius is ~ 6 x 10 8 yr. 
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8 DISCUSSION 

As discussed in Section 4.2, the region of enhanced X-ray 
emission (and enhanced SZ effect; Komatsu et al. 2001) to 
the southeast of the main X-ray peak is significantly hot- 
ter (kT = 18 ± 3 keV) than at the same distance in other 
directions (kT — 12 ± 1 keV). Following Markevitch et al. 
(1999), we can model the properties of this region as a one- 
dimensional shock. In this case, the relative velocity, v, of 
the infalling subcluster and the shock compression factor, 
pi I po, can be estimated from the pre-shock (kTo) and post- 
shock (kTi) temperatures. For kTo = 12 keV and fcTi = 18 
keV, we obtain v — 2250 kms" 1 and pi/ pa = 1.7. If we in- 
stead adopt kTo = 5 keV, a temperature consistent with the 
predicted mass of the merging subcluster from the strong 
lensing analysis discussed in Section 6.2 (using the mass- 
temperature relation of Allen et al. 2001c), then for fcTi = 18 
keV we obtain v — 4550 kms" 1 and pi/ po = 3.0. We note 
that the roughly spherical appearance of the shocked gas 
suggests that this material has expanded following the ini- 
tial shock. Since the material will cool as it expands, it is 
likely that the initial post-shock temperature exceeded 18 
keV and, for kTo = 5, that the relative collision velocity ex- 
ceeded 5000 kms" 1 . Such values are broadly consistent with 
the peak X-ray surface brightness enhancement of a factor 
~ 3 observed at a radius of ~ 180 kpc. 

The second brightest galaxy in the cluster lies approx- 
imately 18 arcsec to the east and 2 arcsec to the south of 
the dominant cluster galaxy. This galaxy has an extended, 
diffuse halo and appears to lie at the centre of a significant 
sub- concentration of galaxies within the cluster (Fig. ^a). 
These properties are consistent with it having been the dom- 
inant galaxy of the recently merged subcluster. The region 
of shocked X-ray gas lies approximately 11 arcsec to the 
south and 6.5 arcsec to the west of the second brightest 
galaxy. Assuming that the dark matter and galaxies as- 
sociated with the subcluster have moved on ahead of the 
shock, it seems plausible that the subcluster travelled in 
from a PA of ~ 210 degrees. This is consistent with the 
extension of the X-ray emission observed on large scales to 
the south (Fig. |l|b) and the weak lensing mass map pre- 
sented by Fischer & Tyson (1997). Given the separation 
of the shocked gas and second brightest galaxy, we esti- 
mate that the shock probably occurred a few 10 7 years ago. 
This timescale is consistent with the spatial extent of the 
shocked gas, assuming that it has expanded at the sound 
speed, c s ~ 1500(T/T0 8 K)" 1/2 kms" 1 . Note also that this 
timescale is comparable to the time required for the elec- 
trons and ions to reach equipartition, so it is possible that 
the electron and ion populations in the shocked gas have 
slightly different temperatures. 

Despite the evidence for shocked gas in the southeast 
quadrant, the Chandra data for the rest of the cluster place 
RXJ1347.5-1145 on the M 25 oo - kT 25 oo and fcT 25 oo - £2500 
relations for relaxed, massive clusters discussed by Allen, 
Schmidt & Fabian (2001). As discussed in Section 4.4, the 
shocked gas produces only a small rise in the mean, emission- 
weighted temperature measured within the central 0.74 Mpc 
radius. The excess bolometric luminosity of the southeast 
quadrant between radii of 60-195 kpc (where the shock is 
most apparent) is ~ 10 45 ergs -1 , which corresponds to only 
~ 5 per cent of the total cluster luminosity (Section 4.4). It 



therefore appears that the overall temperature and luminos- 
ity of the cluster have not been boosted substantially by the 
merger event. This is supported by the consistent X-ray and 
weak lensing mass results discussed in Section 6.1. 

Ricker & Sarazin (2001) and Ritchie & Thomas (2002) 
discuss simulations of mergers between clusters with similar 
central densities and argue that such mergers can lead to 
significant short-term increases in the overall X-ray temper- 
atures and luminosities of clusters during the periods of clos- 
est approach. These effects are most pronounced and long- 
lived in clusters with relatively low central densities (n c < 
few 10" 3 cm" 3 ; Ritchie & Thomas 2002). The fact that large 
boosts in the overall temperature and luminosity are not ob- 
served in RXJ1347. 5-1145 may be related to the very high 
central density in the main cluster core (n e > 10" 1 cm -3 ; 
Fig. |ic|a) and could indicate that the recently- merged sub- 
cluster had a relatively low central gas density and inter- 
acted only weakly with the main cluster core. 

Recently, Cohen & Kneib (2002) have published spec- 
troscopic redshift measurements for 47 cluster members 
for which they determine a velocity dispersion of 910 ± 
130 kms -1 . Inspection of their Fig. 2 suggests that a signif- 
icant fraction of the galaxies with measured redshifts may 
have been associated with the recently merged subcluster. 
This could explain the relatively low velocity dispersion, in 
comparison to the X-ray and lensing results. 

The agreement between the Chandra X-ray and gravita- 
tional lensing mass measurements for RXJ1347. 5-1145 rein- 
forces the results and conclusions drawn from previous stud- 
ies of other relaxed lensing clusters e.g. Abell 2390 (Allen et 
al. 2001b) and Abell 1835 (Schmidt et al. 2001); see also the 
results for MS1358.4+6245 by Arabadjis et al. (2002). The 
close agreement of the independent X-ray and lensing masses 
indicates that the mass measurements are robust and limits 
the contributions from non-thermal sources of pressure sup- 
port in the X-ray gas, such as bulk and/or turbulent motions 
and magnetic fields, in the relaxed regions of the cluster. 



9 CONCLUSIONS 

The main conclusions from this work may be summarized 
as follows: 

(i) We have reported Chandra observations of 
RXJ1347. 5-1145, the most X-ray luminous cluster of galax- 
ies known. We have identified a region of shocked gas (en- 
hanced X-ray brightness and temperature) to the southeast 
of the main X-ray peak, at a position consistent with the 
region of enhanced SZ effect reported by Komatsu et al. 
(2001). The shocked gas probably results from recent sub- 
cluster merger activity. The merger appears not to have 
boosted the overall luminosity and temperature of the clus- 
ter substantially. 

(ii) Excluding the data for the southeast quadrant, we 
have measured the density, temperature and mass profiles 
for the cluster. The mass profile can be parameterized by 
an NFW model with a scale radius r B = 0.40to'.i2 Mpc and 
a concentration parameter, c = 5.9 ± 1.4 (68 per cent con- 
fidence limits). The normalization of the mass profile may 
also be expressed in terms of an effective velocity dispersion, 
a = V50H(z)r s c = 1450^0 kms" 1 . 
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(iii) The best-fit Chandra mass model for RXJ1347.5- 
1145 is in good agreement with independent measurements 
from weak gravitational lensing studies (Fischer & Tyson 
1997). The observed strong lensing configuration in the clus- 
ter core can also be explained with the introduction of an ad- 
ditional mass clump centred on the second brightest galaxy, 
which is likely to have been the dominant galaxy of the re- 
cently merged subcluster. 
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